The Electronic States of Two Oppositely doped Mott Insulators Bilayers 
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We study the effect of Coulomb interaction between two oppositely doped low-dimensional t J 
model systems. We exactly show that, in the one-dimensional case, an arbitrarily weak interaction 
leads to the formation of charge neutral electron-hole pairs. We then use two different mean-field 
theories to address the two-dimensional case, where inter-layer excitons also form and condense. We 
propose that this results in new features which have no analog in single layers, such as the emergence 
of an insulating spin liquid phase. Our simple bilayer model might have relevance to the physics of 
doped Mott insulator interfaces and of the new four layer Ba2Ca3Cu40s compound. 



•4-* 



. 

C ■ 

o ■ 

> ■ 
oo : 

in ■ 
o ■ 

: 

o . 

-i— > ■ 
a ■ 

s; 

T3 ■ 

c : 
o . 
o ■ 



X 



In recent years, the advance in material science made 
it possible to grow atomically smooth epitaxial interfaces 
between different iso-structural transition metal oxides 
[HQ. In particular, it is possible to grow atomically or- 
dered (doped) Mott insulator interfaces 0, 0] ■ Knowing 
the wealth of new physics the interface between (doped) 
band insulators have produced, and the richer physics of 
bulk doped Mott insulators as compared to doped band 
insulators, the interface science of Mott insulators is an- 
ticipated to be a promising direction in years to come. 

Prompted by the above motivation, in this letter we 
take the initiative to explore the interfacial electronic 
states of two low-dimensional doped Mott insulators. We 
concentrate on the specific case where each Mott insula- 
tor is doped with opposite signed charges, and address 
the possible states caused by the Coulomb attraction be- 
tween the resulting n and p type carriers. In particular, 
we find that the formation of electron-hole pairs, also 
known as excitons, across the interface gives rise to new 
physics which is absent in the isolated systems. 

More specifically, below we study the interface between 
either two one-dimensional (ID) tJ model chains or two 
two-dimensional (2D) t J model layers, where the upper 
chain (or layer) is doped with a density x of electrons and 
the lower chain (or layer) has the same density of holes. 
We first consider the ID case, where we explicitly show 
that an arbitrarily small inter-chain Coulomb interaction 
leads to exciton formation. We then consider the exper- 
imentally relevant 2D system which, aside from being 
an explicit model for engineered (doped) Mott insulator 
interfaces 0, 0|, may also be pertinent to address the 
recently studied stoichiometric four-layer Ba 2 Ca3Cu40s 
compound where the simultaneous presence of n and p 
type Fermi surfaces indicates the existence of a "self- 
doping" effect 0. Unfortunately, the exact argument 
used in one dimension does not apply in higher dimen- 
sions and, in this case, we rather resort two different 
mean-field schemes, namely, a recently developed fully 
fcrmionic mean-field theory [5L JfJ and the traditional 
slave-boson mean- field theory [7], [8( . Although these the- 



ories sometimes miss the ultimate low energy physics, 
they have successfully predicted the qualitative nature of 
the cuprate phase diagram [§[ • Since our current pur- 
pose is to explore the possible interface electronic phases, 
we feel that it is meaningful to work out the prediction 
of the aforementioned approaches, which we find yield 
consistent results. 

The model. The Hamiltonian describing our interface 

is Hi n terface = Yl<a=u,d + H int where 



(ij)eNN 



) - Kc\,i<3 + h.c) (1) 



is the Hamiltonian for the upper (a = u) and lower 
(a = I) systems (chains and layers in the ID and 2D 
case, respectively). J and t are the nearest neighbor 
(NN) antiferromagnetic exchange and electron hopping 

J 



parameters, = V[c\ i V c^JP, where c' a i (y creates 
electrons in chain (layer) a. V projects out doubly occu- 
pied (vacant) sites when a = I (a = u). n a ^ = c' a ^. a ^ and 

S a ,i — hc a i<yc a ,% (<J are the Pauli matrices) are the elec- 
tron number and spin operators. The across the interface 
interaction is described by 

Afj, 



H, 



V}(l-n u<i )(l-n lti ) + 



where V > is the Coulomb repulsion parameter and 
Afi is the energy cost to move an electron from the lower 
to the upper chain (layer). 

The exactly solvable ID interface. We first look at 
the interface model for two ID chains. The model then 
has the remarkable property that in the limit J — > 
the wavefunction of the ground states factorizes into a 
spin wavefunction and a charge wavefunction. This in- 
teresting feature was first observed in Ref. [l(| for the 
integrable Hubbard chain at large U. However, it is also 
common to various non-integrable t J type models, lead- 
ing to a variety of exact statements [ll|, [l2| . We proceed 
by performing a particle-hole transformation on the up- 
per chain, which makes doped electrons formally iden- 
tical to holes, while changing the sign of the interchain 
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interaction, V — > — V. We then introduce the wavefunc- 
tion tp(xi,ai . . .x N ,a N ; yi,xi ■ ■■Vn,Xn), where x u a l 
are the positions and spins of the particles in the upper 
chain, and y^, \i those of the lower chain. Both sets of co- 
ordinates Xi and yi are assumed to be given in ascending 
order. The action of the Hamiltonian on the wavcfunc- 
tion ip now takes the following form 

H^= -t^ =1 ^(...xi±lai...)+^(... m ±l X i...)) 



-VEL EZ =1 8*„y m <!>(■■■) + Hj H- ■■) (3) 

where Hj represents the spin exchange term proportional 
to J. In the limit J — > one can easily verify that the 
eigenstates of the system described by Eq. (j3|) are of the 
following factorized form: 

. . ) = f(xi ...x N ]y x ... y N ) g(ai ...a N ;xi--- Xn)- 

(4) 

This follows from the fact that at J = 0, the Hamilto- 
nian Eq. |3]) docs not act on the spin degrees of free- 
dom. Hence, / must be an eigenfunction of Eq. (|3|) with 
J = 0, whereas g can be an arbitrary function of spin 
variables. Precisely at J = the ground state is thus 
hugely degenerate in the spin sector. This degeneracy 
can be lifted in first order perturbation theory via an ef- 
fective spin Hamiltonian acting on a spin ladder of 2N 
spins, where N is the number of singly occupied sites in 
both the upper and the lower chain, and the holes are 
effectively "squeezed" out of the system [ll], [ijj . At first 
order in J, the spin wavefunction g is thus uniquely deter- 
mined. Here we are, however, more concerned with the 
charge correlations of the system. In this regard, we ob- 
serve that the determination of the charge wavefunction 
/ from the Hamiltonian Eq. ^ is equivalent to finding 
the ground state of the attractive Hubbard model. In- 
deed, the effective charge Hamiltonian acting on / can 
be rewritten in second quantized form as 

H c = -t^(a\ +l>a a ita +h.c)-V^a\ tU ai tU a\ d ai4 (5) 

where the fermion operator ai^ a carries a "pseudospin" 
index a labeling the layer. It is well known that Eq. (JSJ) 
is in the spin gapped Luther-Emery liquid phase of Id 
systems for any V > 0, which is the ID analog of a paired 
superfluid [ll[ . The spin gap of the Hubbard model hence 
implies a charge gap due to the formation of electron- 
hole pairs in the double-chain system. The important 
lesson to be learned is that, in one dimension, the phase 
diagram of our interface model features exciton formation 
at arbitrarily weak inter-chain interaction. 

The fully fermionic mean-field theory. We now 
discuss the 2D bilayer Hamiltonian Hi nte rface at the 
mean-field level. We start with the fully fermionic mean- 
field theory introduced in Refs. 0, @. This theory in- 
troduces fermionic operators d\ i = [d\ i ^d\ i J that cre- 
ate charge (+e) and spin- 1/2 holes in the lower Hubbard 
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band of the lower layer (a = I) and charge (— e) and 
spin-1/2 electrons in the upper Hubbard band of the up- 
per layer (a = u). The Mott insulating background on 
top of which these carriers move is described by a singly 
occupied lattice of neutral spin-1/2 fermions whose cre- 
ation operators are = [f^ i + f a i J. In terms of the d 
and / fermions the projected electron operators can be 
expressed as Q 
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9 s " d 1 
ti,i u,i,-(T 



(6) 



where s a = (+1), (-1) for a =t,4, S a ,i = jfl^fa.i and 
d a ,i,a = <4,i,o-(l - d a ,i,-crda,i,-o-)- These relations allow 
us to recast Hi nter face in terms of the d, d) and /, f' op- 
erators, after what we perform a mean-field factorization 
to obtain a quadratic Hamiltonian. 

First, we summarize the (single layer) mean- field re- 
sults of Ref . [H, [|| . At low enough doping the t J model 
develops antiferromagnetic order [13j . which yields non- 
zero staggered magnetization order parameters 



i/=2,3 



-1) 

■ E 

u£vNN 



vhva,i+u + h.c.) (7) 



Here, ^ = [< M <, 2 ] = s^[fl, A fa^i] and ^ = 

[vl,i,ivl,i,2\ = 8 a~ %v [di M d atill ]; s a = (+1),(-1) for 
a = l,u; ti = J and tj, — J/2: u — ±x ± y for v = 2 
and u = ±2i, ±2y for v = 3 [141 ] . In addition to the an- 
tiferromagnetic order parameter, singlet bond order pa- 
rameters which capture the spin liquid correlations are 
also introduced 

x = ^U^a tj ) ; A = (-) W, (<^ aj ) (8) 

with (ij) S NN. Finally, the Bose amplitudes which de- 
scribe the quantum coherence between the doped carriers 
and the spin background are given by 



E 

w=l,2,3 u£vNN 



(9) 



where t\ = t and u = ±x, ±y for v = 1. 

In terms of these order parameters, the in-plane mean- 
field Hamiltonian is given by [B|, Q 



K F = E [ < k < 

a,k 



E 

a . k 



a z k a z + al<r x (i k a z 



^a,k + ^IkVlk+QVa.k 



tpa,k 
Va,k 



(10) 



In the above equation a z k 



(¥x - hx)(coak x + 



cos ky ) + ao where qq is the Langrange multiplier that 
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FIG. 1: tj model bilayer mean-field phase diagram for t — 
3 J. It includes the antiferromagnetic insulator (AFI), the 
spin liquid insulator (SLI), the superconductor (SC), the band 
insulator (BI) and the Fermi liquid metal (FLM) described in 
the main text. 



ensures {f ai f a ,i) = 1, and J = (1 — x) 2 J. In addi- 
tion, a% = — ^fA(cosk x — cos k y ), Pk = — ^g 2 - [h (cos k x + 
cosky) + 2ti cos k x cos k y + t 3 (cos2k x + cos2fc y )] + b\, 
7fc = h cos k x cos k y + Y( cos ^k x + cos2fc y ) + vf = 
2 Sa (n - 0.34Jto), ^ = 2s a m(^ - T £). , 

Unlike other slave-particle approaches [7], la |l5j, all 
particles in the above mean-field theory are fermions. 
The central advantage of this mean-field approach is the 
ease to describe fermionic quasiparticles, which are cre- 
ated by the eft operator. The "holon" in the slave-boson 
approach, which is a spin singlet and charge +e non- 
quasiparticle object, is the composite bosonic particle 
made up of a d quasiparticle and a / neutral fermion, and 
its condensation is described by the the order parameters 
fro,i @- The fully fermionic approach, which yields re- 
sults consistent with numerical calculations Q, success- 
fully reproduces various non-trivial aspects of ARPES 
and STM experiments in the cuprate superconductors, 
including the Fermi arcs, the dispersion kink, the min- 
imum gap locus, the flat dispersion around (7T,0), the 
peak-dip-hump and the absence of coherence peaks in 
deeply underdoped dl/dV spectra 0, Encouraged 
by the above successes we extend Eq. (ITU1) to describe 
two oppositely doped Mott insulator layers. Following 
our exact ID result, we introduce the order parameters 

e * = (vliV-Vuj + vl,i<r-Vi,i) 

e* = (-)^ +J » (vli<r-Vu,i - vlj<r-m,i)> ( n ) 

to describe the formation of spin singlet and spin triplet 
excitons. In Eq. (fTT|l er_ = (a x — ia y )/2. The resulting 
inter-layer mean-field Hamiltonian reads 

H^ t F = — o £ ( e sVl, k <rxVl,k + ietvl k+Q VyVLk + h.c)j 

k 

(12) 

The zero temperature mean-field phase diagram is de- 
termined upon minimizing the mean-field energy (H^ IF + 

Hf F + Ht, 1 / - E ), where E = -[Hf{x 2 + A 2 ) + 



IMJNm 2 - 8Nmn - ANboh + ^N{e% + e|) + A//JV] 
and N is the number of lattice sites. Motivated by the 
mutual exclusion of antiferromagnetism and supercon- 
ductivity [H in the cuprates, we omit the possible coex- 
istence of antiferromagnetic and superconducting orders. 
The resulting phase diagram for the bilayer t = 3 J model 
is depicted in Fig. [TJ It comprises five distinct regions 
whose properties we discuss below: 

(i) AF insulator (x,m,n,e s ,e t ^ and &o,foi — 0). 
Since m ^ spins are antiferromagnetically ordered. 
Electrically neutral excitons form (e s , ej ^ 0) which cost 
a finite energy to break, thus opening a gap in the charge 
sector. 

(ii) Gapless spin liquid insulator (%, A, e s ^ and 
to, n, et,bo,b\ = 0). In the absence of AF order [m = 0), 
the spins are in a gapless d-wave (A ^ 0) paired spin liq- 
uid state reminiscent of the pseudogap state in the single 
layer problem ||. However, in the bilayer case excitons 
(e s 7^ 0) render the system a charge insulator. 

(Hi) (D + S)-wave superconductor (x, A, bo, &i, e s ^ 
and rn, n, e% — 0). Doped carriers hybridize with the d- 
wave paired spins (A, bo ^ 0) and, at the same time, 
participate in forming s-wave excitons. As a result, the 
superconducting gap acquires a D + S symmetry. 

(iv) Band insulator [\, bo, b\, e s ^ and to, n, et, A = 
0). In the absence of pairing (A = 0) there is no super- 
conducting order. The non-zero &o and e s condensates 
open a gap in both the spin and charge sectors and the 
system behaves as a renormalized band insulator. 

(v) Fermi liquid metal (x,foo ; ^i ^ an d 
m,n,et,e s , A = 0). In the absence of A and e s 
the quasiparticle gap closes and the system behaves as a 
renormalized Fermi liquid. 

The phase diagram in Fig. [JJ results from the competi- 
tion between the spin exchange energy, the kinetic energy 
(of the doped carriers) and the inter-layer Coulomb en- 
ergy. The latter two vanish at half-filling, hence, at suffi- 
ciently low doping, antiferromagnetic order prevails. As 
doping increases, the itinerancy of charge carriers tends 
to destroy antiferromagnetism and to favor a spin liquid 
state instead [ljj • In that case, and for sufficiently large 
V , the antiferromagnetic insulator is replaced by the spin 
liquid insulator where all doped carriers form inter-layer 
excitons. Since charge carriers are more mobile when 
they are not part of an exciton, partial exciton unbinding 
is favored at low V . The carriers thus liberated from the 
excitonic condensate then bind to the spins and, in the 
presence of spin pairing, the system is a superconductor. 
At large doping, kinetic energy dominates over the ex- 
change interaction and, as a result, the superconducting 
gap closes, rendering the bilayer metallic. This metallic 
state occurs in a system with two electrons per bilayer 
unit cell and, nominally, it is rather expected to be an 
insulator. Such band theory expectations are amiss due 
to the effect of strong correlations [l8j]. Since the effect 
of electron-electron correlations decreases upon further 
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doping, the system eventually becomes a band insulator, 
as expected in the weak coupling limit. 

The slave-boson mean-field theory. To vindicate 
the above predictions we also determine the phase dia- 
gram using the slave-boson mean-field theory 0, [|| . We 
use the antiferromagnetic (to), in-plane d-wave pairing 
(A-f), inter-plane exciton pairing (A h ) and single boson 
condensation ((&)) mean-field order parameters and, as 
before, ignore the coexistence of antiferromagnetic and 
superconducting orders. To our satisfaction, the results 
are in qualitative agreement with those obtained using 
the above fully fermionic mean-field theory. Upon in- 
creasing the carrier concentration, the excitonic antifer- 
romagnetic insulator (to 0, A b ^ 0, A f ^ 0, (6) = 0) 
present at < x < x c evolves to the excitonic d-wave 
superconductor (m = 0, A b ^ 0, A* ^ 0, (b) ^ 0) for 
x > x c . The critical concentration x c increases with V/J 
in a manner consistent with the AFI-SC phase bound- 
ary depicted in Fig. [TJ If V/J is larger than a certain 
value, however, we find that the phase that replaces the 
AF state is a spin liquid insulator (to = 0, ^ 0, 
A 6 7^ 0, (6) =0). A similar AFI-to-SLI transition is 
found for large V/J values in Fig. [TJ Detailed findings of 
the slave-boson calculation are presented elsewhere [19j ■ 

Discussion. Above, we use exact arguments and two 
different mean-field theories to conclude that Coulomb 
interaction between oppositely doped tJ model systems 
leads to the formation of excitons in ID and 2D, respec- 
tively. Since in the ID [2D] system the hole momentum 
distribution is shifted from that of electrons by tt [(tt, it)], 
excitons carry a finite momentum equal to tt [(tt, tt)] [141 ]. 

Even though in the 2D case our model Hamiltonian 
Hinterface is, at most, a simplified version of material 
compounds' Hamiltonians, it provides the playground to 
explore the important effect of inter-layer Coulomb in- 
teraction, which is the largest interaction between neigh- 
boring oppositely doped Mott insulators. We also believe 
that the single layer physics is correctly captured by re- 
stricting the in-plane Coulomb interaction to electrons 
sharing the same lattice site. However, it is well known 
that terms describing electron hopping beyond NN sites 
can significantly alter the single plane properties 17 , 2(| ■ 
Above, we neglect such terms, which may limit our abil- 
ity to address systems like the P^CasCruOs compound 
where farther neighbor hopping processes are believed to 
be important. Although inclusion of longer range hop- 
ping is the subject of future work, we expect the present 
conclusions to apply as long as xt/t', xV/t' > 1. 

Finally, we note that Hi nter f ace is interesting in itself, 
as it illustrates how exciting new physics (absent in iso- 
lated single layers) may emerge in bilayer systems. In the 
present case, such physics is summarized as follows: 

(i) In the presence of large inter-layer Coulomb interac- 



tion all doped carriers form electrically neutral electron- 
hole pairs and the system is an insulator. As these pairs 
move throughout the layers, they disorder the spin back- 
ground and stabilize a fractionalized spin liquid state 
which supports neutral spin-1/2 excitations. So far, there 
is no evidence for such spin liquid states in single layers 
2l|, where the itinerant charge carriers render the spin 
liquid unstable to superconductivity at low temperature. 
Here, we propose that the charge gap in the excitonic in- 
sulating phase can protect the exotic spin liquid state in 
the bilayer formed by oppositely doped Mott insulators. 

(ii) The coexistence of the excitonic condensate and su- 
perconductivity changes the quasiparticle gap from pure 
D to D+S symmetry. The extra gap component deviates 
the nodes away from the (0, 0) — (ir, tt) and (0, 0) — (— tt, tt) 
directions. This effect can, in principle, be detected in 
ARPES experiments. 
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